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^Bottom-loss  measurements  were  made  for  frequencies  less  than  800  II:  at  three 
stations  in  the  llattoras  .\byssal  Plain.  The  stations  were  each  separated  approx- 
imately 55  km  alons  70®30’W  longitude.  The  data  acquisition  procedure  was  to 
launch  two  sol f-recordins  submersibles  (Al/TOBUOYS)  and  then  maneuver  from  aiul  to 
the  launch  site  detonatins  explosives  as  the  AUTOBUOYS  recorded  data.  Because  ot 
the  experimental  geometry,  the  direct  and  bottom-interacting  signals  wore  time- 
separated  as  they  arrived  at  the  AUTOBUOYS.  Thus,  a "self-calibrating"  (^over) 


20.  (Cont'd) 

i processing  procedure  could  be  used  to  compute  bottom  loss  values. 

Bottom  loss  versus  grazing  angles  (5-66  deg)  was  obtained  over  a wide 
bandwidth  (90-790  Hz)  and  for  a number  of  narrow  bandwidths  (1/3  octave  and 
less).  Broadband  results  were  found  to  be  homologous  between  stations, 
whereas  the  narrowband  results  were  less  correlated.  Negative  bottom  loss 
values  were  also  observed. 

Replica  deconvolution  processing  was  performed  using  the  associated 
time-separated  direct  and  bottom  interacted  signals  over  a large  range  of 
grazing  angles.  A history  of  the  bottom/subbottoro  Impulse  response  as  a 
function  of  grazing  angle  was  then  constructed  from  the  deconvolved  signals. 
This  history  showed  evidence  of  strong  subbottom  reflected  and  sediment 
refracted  arrivals  in  the  bottom  interacting  signals.  Since  these  types  of 
arrivals  were  present,  the  calculation  of  negative  bottom  loss  can  be 
attributed  to  use  of  the  Rayleigh  plane  wave  reflection  model,  which  does 
not  account  for  non-plane  wave  refracted  or  reflected  signals. 

Analyses  of  the  impulse  responses  for  the  southern  station  also  aided 
in  developing  a geophysical  description  of  the  sediment  structure.  The 
derived  structure  consists  of  a thin  (21  m),^ constant  sound-speed  layer 
overlaying  a thick  layer  having  a linear  depth-dependent  sound  speed.  This 
structure  is  representative  of  a class  of  documented  subbottom  sediment 
structures  that  can  support  the  propagation  of  low-frequency  acoustic 
signals  traversing  non-plane  wave  reflected  and  refracted  acoustic  paths, ^ 

i \ 

\ 


TR  5781 


TABLE  OF  CONTENTS 

Page 

LIST  OF  ILLUSTRATIONS ii 

LIST  OF  TABLES iii 

INTRODUCTION 1 

BASIC  OBJECTIVES.  1 

BASIC  APPROACH 1 

DATA  ACQUISITION  PROCEDURES  2 

Bottom  Loss  Measurements  With  AUTOBOUYS.  2 

Seismic  Measurements  3 

Sediment  Core  Measurements  , 3 

SUMMARY  OF  OPERATIONS  3 

BOTTOM  LOSS  DATA  PROCESSING  PROCEDURES 4 

Digitization  .....  4 

Preprocessing  Digital  Filtering 4 

Processing  for  Separated  Arrivals.  4 

BOTTOM  LOSS  RESULTS  8 

Nature  of  the  Bottom  Arrivals 8 

Quality  of  Data 12 

Geometry  Dependence  of  Bottom  Loss  14 

Frequency  Variations  of  Bottom  Loss 14 

Negative  Bottom  Loss 14 

Geographical  Variations  of  Bottom  Loss < . ..  15 

PRELIMINARY  MODELING  AND  INTERPRETATION  15 

Empirical  Modeling  Via  Polynomials ' . 15 

Input  Parameters  for  a Geophysical  Model  16 

SUMMARY  AND  CONCLUSIONS  17 

REFERENCES 34 

APPENDIX  A--SOUND  SPEED  VERSUS  DEPTH  PROFILES  A-1 

APPENDIX  B— BOTTOM  LOSS  VERSUS  GRAZING  ANGLE  DATA 

FOR  ALL  STATIONS  AND  FREQUENCIES B-1 

APPENDIX  C— SIXTH-ORDER  POLYNOMIAL  REPRESENTATIONS  OF 

BOTTOM  LOSS  DATA C-1 


TR  5781 


LIST  OF  ILLUSTRATIONS 


Figur**  Page 

1 Operations  Area  for  Hatteraa  Abyssal  Plain  Low 

Frequency  Bottom  Loss  Measurementa  .....  19 

2 Typical  Repreaentat ion  of  Geometry  for  Bottom  Loss 

Measurements  (Station  BI.-l)  20 

\ Block  Diagram  of  Data  Processing  System  for  Bottom  Loss 

Measurements  ...........  21 

A Sequence  of  Arrivals  for  a Typical  Reception  Before 

and  After  Filtering  Line  Components  - (Station  BL- I , 

AinXlBlKlV  No.  2,  0j^  - ll.O*') 22 

5 Averaged  Direct  Path  Knergy  Spectrum  (Average  Source 

Depth  909  m)  . . . . .* 2J 

b Example  of  Signal  Enliancement  Obtained  Using 

Deconvolution 2^ 

7 Consecutive  Time  Displays  of  Bottom  Arrivals  as  a 

Function  of  Gracing  Angle  With  Respect  to  Interface 
Reflected  Angle  (Sta^iol^  BL~l) 25 

8 Squared  Travel  Time  Difference  Between  Subbuttom  and 

Ocean  Bottom  Reflections  as  a Function  of  the  Cosine 
Squared  of  the  Ocean  Bottom  Gracing  Angle  2b 

9 Measvired  and  Predicted  Travel  Time  Differences  Between 

Ocean  Bottom  Reflected  and  Subbottom  Refracted 

Arrivals 27 

10  Deconvolved  Sediment  Interacted  Signals  Low 

Pass  Filtered  at  Three  Frequencies;  500,  250, 

and  200  He  (Site  BL- 1 , Upper  Buoy) 28 

11  Comparison  of  Bottom  Loss  Versus  Gracing  Angle  as  a 

Fur  ' a of  Receiver  Depth  (Station  BL-1) 29 

12  Frequency  Variation  of  Bottom  Loss  Versus  Grazing 

Angle 30 

13  Geographical  Variations  of  Bottom  Loss  Versus 

Grating  Angle  Data  (Broadband)  ....  31 


11 


TR  5781 


LIST  OF  ILLUSTRATIONS  (Covtl'd) 


Page 


Polynomial  Fitted  and  Averaged  Broadband  Bottom 
Loss  Meaaured  at  Three  Stations  in  the 

Hatteras  Abyssal  Plain  32 

Geographical  Variations  in  Bottom  Loss  (Narrowband)  33 


LIST  OF  TABLES 


Page 


Locations  and  Descriptions  of  Bottom  Loss 

Measurements 2 

Summary  Analysis  of  Sediments  Recovered  by 

Boomerang  Core  Near  Station  BL-3 !) 

Summary  Analysis  of  Sediments  Recovered  by 

Boomerang  Core  Near  Station  BL-2 5 

Filter  Parameters  Used  In  Narrowband  Bottom  Loss 

Analysis  7 


iii/lv 
Reverse  Blank 


J.  l-  , _ 


TR  5781 


HATTERAS  ABYSSAL  PLAIN  LOW-FREQUENCY 
BOTTOM  LOSS  MEASUREMENTS 


INTRODUCTION 


In  May  and  June  1973,  low-frequency  bottom-loss  measurements  were  con- 
ducted at  three  locations  in  the  Hatteras  Abyssal  Plain  of  the  Atlantic 
Ocean  by  the  Acoustic  and  Environmental  Research  Division  of  the  Naval 
Underwater  Systems  Center  (NUSC)  and  the  Lament -Doherty  Geological  Observa- 
tory (LDGO).  This  project  was  designated  MAINLOBE  (MAjor  INvestigation  for 
Low-Frequency  Ocean  Bottom-Loss  Experiment)  73.  The  text  that  follows  is 
divided  into  the  topical  areas  of  (1)  basic  objectives,  (2)  basic  approach, 
(3)  data  aquisition  procedures,  (A)  summary  of  operations,  (5)  data  process- 
ing procedures,  (6)  analysis  and  interpretation  of  the  bottom-loss  measure- 
ments, (7)  preliminary  modeling,  and  (8)  summary  of  results. 

BASIC  OBJECTIVES 

The  primary  objectives  of  the  experiments  described  in  this  report  arc 
the  (1)  measurement  of  ocean  bottom  loss  over  a wide  frequency  range 
(20-1000  Hz)  as  a function  of  grazing  angle  (5-90  deg)  and  (2)  estimation  of 
sediment  acoustic  parameters  (such  as  density,  sound  speed,  and  attenuation) 
for  sediment  interaction  modeling  and  prediction  of  propagation  loss. 

BASIC  APPROACH 

Three  measurement  locations,  shown  within  the  circled  areas  in  figure  1* 
as  stations  BL-1,  -2,  and  -J,  were  employed  in  MAINLOBE  to  account  for 
possible  acoustic  variations  from  changes  in  bottom  topography  and  structure 
(both  subbottora  and  bottom  interface).  Data  were  collected  successfully  at 
each  of  the  three  measurement  locations.  Explosive  signals,  underwater 
sound  (SUS)  charges  were  used  in  conjunction  with  two  self-contained,  com- 
pletely automated  receiving  systems  (AUTOBUOYS)  to  generate  and  receive  the 
low-frequency  wide-band  underwater  signals.^ 

To  enable  estimation  of  sediment  acoustic  parameters,  seismic  reflec- 
tivity measurements  and  subbottom  profiling  were  planned  along  the  tracks 
shown  as  dashed  lines  in  figure  1.  These  operations,  conducted  jointly  by 
personnel  from  NUSC  and  LDGO,  involved  the  transmission  of  acoustic  signals 
generated  by  air  guns  and  the  reception  and  recording  of  these  signals  by 
means  of  a towed  horizontal  receiving  array  system. ^ Boomerang  cores  were 
also  planned  at  locations  near  the  bottom-loss  measurements.  The  positions 
of  the  two  successful  boomerang  cores  are  shown  in  figure  1. 

♦Figures  1 through  15  are  presented  at  the  end  of  the  text. 
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Envi ronaental  and  oceanographic  data  in  aupport  of  the  above  aeaaure- 
menta  included: 

1.  found  apeed  profilea  (SSP)  (atationa  BL~li~2,~3)  and 

2.  expendable  bathytheraographa  (XBT) 


DATA  ACQUISITION  PROCEDURES 


BOTTOM-LOSS  MEASUREMENTS  WITH  AUTOBUOYS 

The  ocean  bottom-loaa  experiuenta  conducted  during  project  MAINLOBE 
involved  the  deployment  of  two  AUTOBUOYS,  SUS  exploaive  chargea,  and 
AM/SSQ-37A  or  AN/SSQ-41  sonobuoya.  An  example  cf  the  meaauremcnc  geometry 
is  shown  in  figure  2 for  a typical  AUTOBUOY  receiver.  A primary  objective 
in  NUSC'a  approach  to  measurement  of  low-frequency  acoustic  bottom  loss 
using  SUS  charges  is  to  obtain  t ime- isolated  direct  and  bottom  interaction 
path  arrivals  over  a wide  range  of  graxing  angle.  The  optimum  depths  of 
signal  source  and  receiver  for  arrival  separation  were  determined  from 
ray-tracing  analysis  using  the  Navy  Interim  Surface  Ship  Model  II  (NISSM  II) 
computer  model. ^ To  study  the  possible  effects  of  source/receiver 
geometry  on  bottom  loss  measurements,  two  AUTOBUOYS  were  deployed  at  each 
measurement  station.  The  general  measurement  procedure  consisted  of  (1) 
deploying  the  first  AUTOBUOY  on  station,  (2)  traversing  2.3  km  north,  (3) 
deploying  the  second  AUTOBUOY  on  station,  (4)  traversing  32  km  north,  (5) 
turning  180  deg,  and  (6)  deploying  SUS  charges  while  returning  on  a due 
south  track  to  the  second  AUTOBUOY.  The  SUS  charges  were  type  EX-94(0), 
containing  816  g tetryl  and  31  g primer,  set  to  detonate  at  914  m. 

Data  concerning  AUTOBUOY  hovering  depths,  functioning,  and  minimum/ 
maximum  range  separation  are  listed  in  tablo  1 for  all  three  measurement 
stations.  The  sonobuoys  were  used  to  monitor  the  source  level  and  detona- 
tion depth  of  the  explosives  and  the  near  range  (high  grazing  angle)  bottom 
path  arrivals  during  the  measurement  period. 

Table  1.  Locations  and  Descriptions  of  Bottom  Loss  Measurements 
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E.<i  leceivet  depths  see  Appendti  A. 


Either  immediately  prior  to  or  following  the  acoustic  measurements,  a 
deep  SSP  cast  was  made  for  each  station  to  a depth  of  at  least  4300  m. 
These  sound  speed  versus  depth  profiles  are  presented  in  appendix  A. 
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SEISMIC  MEASUREMENTS 

The  procedure  for  obtaining  wide  angle  aeiamic  data  wae  to  first  deploy 
an  air  gun  and  then  pay  out  a buoyant  aeitmic  receiving  line  array  (122  m 
long)  from  the  stern  of  the  research  ship  as  it  proceeded  at  4 knots.  After 
deployment,  the  (peed  of  the  ship  was  increased  to  6 knots  and  a sonobuoy 
was  deployed.  The  acoustic  signals  generated  by  the  air  gun  were  received 
by  the  line  array  and  the  sonobuoy  and  recorded  on  LDCO's  seismic  profiling 
system.  The  profiling  system  provides  (in  situ)  some  general  results  con- 
cerning the  characteristics  of  the  ocean  sediment  and  layer  structure.  The 
same  data  were  also  recorded  on  analog  tape  for  further  laboratory  analy- 
sis. Additional  information  concerning  the  LDCO  seismic  data  processing 
procedures  can  be  found  in  "Physics  of  Sound  in  Marine  Sediments."^  An 
analysis  of  measurement  results  is  reported  in  an  unpublished  work  by  Houtz 
and  Bryan. ^ 

A buoy  having  a deep  (1128  m)  hydrophone  was  also  deployed  while  making 
the  seismic  measurements.  This  acoustic  radio  buoy  (ARAB)  is  a NUSC  innova- 
tion consisting  of  sonobuoy  transmitters  installed  in  a spar  buoy  housing. 
The  ARAB  system  was  used  primarily  to  obtain  a measure  of  the  effective 
low-frequency  beam  pattern  (including  the  Lloyd  mirror  effect)  of  the  air 
gun.  Poor  signal-to-noise  ratio  (SNR)  of  the  bottom-interacted  signals 
precluded  analysis  of  the  ARAB  data  to  obtain  high  grazing  angle  bottom  loss 
resul ts . 


SEDIMENT  CORE  MEASUREMENTS 

Core  measurements  were  made  using  a Benthos  boomerang  corc.^  This 
device  is  an  untethered  free-fall  corer  having  a 1.22  m long  barrel  and  an 
outside  diameter  of  7.94  cm.  The  plastic  core  liner  tube,  which  holds  the 
core  sample,  is  1.22  m long,  has  an  outside  diameter  of  7.3  cm,  and  an 
inside  diameter  of  6.50  cm.  After  the  core  penetrates  the  sediment,  a 
trigger  releases  the  barrel  from  the  housing  and  buoyant  glass  spheres  lift 
it  to  the  surface  with  the  captured  sediment. 

SUMMARY  OF  OPERATIONS 


Only  a summary  of  project  MAINLOBE  operations  is  presented  herein 
because  daily  operations  have  been  reported  in  detail.^  Successful 
bottom-loss  measurements  were  made  at  all  three  stations:  however,  no 
usable  data  were  obtained  from  AUTOBUOY  No.  1 at  stations  BL-2  and  BL-3 
because  it  malfunctioned. 

Seismic  measurements  were  successfully  completed  along  the  tracks  shown 
in  figure  1.  The  ARAB  was  successfully  deployed  in  each  of  the  three  areas 
of  seismic  profiling;  however,  low-frequency  noise  caused  by  electrical 
feedback  in  the  ship  power  line  during  operation  of  the  air  gun  compressor 
degraded  the  data.  Seismic  data  were  not  affected  by  this  problem. 
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Only  two  of  six  boomerang  cores  were  successfully  retrieved.  The  first 
wus  retrieved  at  29®  06.8'N,  70®  30.0'W,  near  station  BL-3.  The  other 
was  retrieved  at  28®  36.5'N,  70®  24.5'W,  near  station  BL-2.  A summary 
analysis  of  the  data  is  presented  in  tables  2 and  3. 

P'.n'OM-LOSS  DATA  PROCESSING  PROCEDURES 

The  basic  approach  to  the  measurement  and  processing  techniques  employed 
at  NUSC  for  the  determination  of  ocean  bottom  loss  has  been  previously  docu- 
mented. A suimnary  of  these  procedures  and  the  special  processing 

techniques  required  is  described  below. 


DIGITIZATION 

As  shown  in  figure  3,  data  from  two  of  the  three  available  staggered- 
gain  AUTOBUOY  analog  tape  channels  were  passed  through  a 750-Hz  low-pass 
filter,  digitally  sampled,  and  stored  on  magnetic  tape  using  a UNIVAC  1230 
computer  system.  To  compensate  for  record  and/or  reproduce  tape-speed 
variations,  external  sampling  commands  were  generated  from  a 1000-He  time 
code  carrier  (also  recorded  on  the  AUTOBUOY  analog  data  tape)  via  a phase- 
locked  frequency  multiplier.  A sampling  rate  of  3000  Hz,  which  exceeds  the 
Nyquist  sampling  criterion  rate  by  a factor  of  2,  was  used  to  facilitate 
deconvolution  processing  of  the  data. 


PREPROCESSING  DIGITAL  FILTERING 

Following  digitization,  contaminating  interference  caused  by  alternating 
current  components  during  playback  of  the  AUTOBUOY  analog  tape  was  removed 
by  digital  filtering.  A digital  filter  was  designed  to  pass  all  frequencies 
except  the  narrowband  alternating  current  line  components.  In  figure  4,  a 
typical  sequence  of  arrivals  prior  to  and  following  the  line  component 
filtering  is  shown. 


PROCESSING  FOR  SEPARATED  ARRIVALS 


Bottom  Loss  Calculation 

The  determination  of  bottom  loss  for  any  given  explosive  detonation  was 
accomplished  by  using 

®L  “ ^*^D  “ *^B^  * (*^D"^B^» 

where 

B^  <■  bottom  loss  (dB), 

Rd  ■ received  energy  in  direct  arrival  (dB), 
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Table  2.  Siiinmai'y  Ai.alysis  ov  Sediments  Recovered 
by  boomerang  Cere  Near  Station  BL-3 
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Depth 
( cm) 

Sound  Speeds 
(m/  s ) 

Bulk  Density 
(g/cm^) 

Impedance 
( 10^. g/ cn^/ s) 

Inter  face 

1608 

1.11 

1.785 

15.2 

1600 

1.11 

1.776 

30.5 

1637 

1.11 

1.817 

45.7 

1603 

1.15 

1.843 

57.2 

1637 

1.15 

1.883 

Average 

1617 

1.13 

1.827 

♦Corrected 

to  seafloor  in  situ 

conditions. 

Table  3.  Summary  Analysis  of  Sediments  Recovered 
by  Boomerang  Core  Near  Station  BL-2 

Depth 
( cm) 

Sound  Speed* 

(ra/  s) 

Bulk  Density 
(g/ cm^) 

Impedance 
(10^  g/cm^/s) 

Inter  face 

1609 

1.33 

2.140 

3.8 

1584 

1.33 

2.107 

12.7 

1616 

1.33 

2.149 

29.4 

1616 

1.33 

2.149 

41.2 

1621 

1.20 

1.945 

Average 

1609 

1.30 

2.092 

♦Corrected 

to  ‘^eafloor  in  situ 

conditions . 
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Rg  » received  energy  in  bottom- interacting  arrival  (dB), 

Np  ■ propagation  loss  over  direct  path  (dB),  and 

Ng  •"  propagation  loss  over  bottom-reflected  path  (dB). 

As  is  evident  from  the  above  expression,  the  procedure  for  calculating 
bottom  loss  is  comparative  and  self-calibrating  in  that  only  the  relative 
energy  difference  between  the  direct  and  bottom  interacting  arrivals  need  be 
experimentally  determined.  The  corresponding  propagation  loss  within  the 
water  column  for  the  two  arrivals  is  determined  using  a NISSM  II  computer 
ray-tracing  program.^  The  input  parameters  for  NISSM  II  were  (1)  in-situ 
measured  water  SSP,  (2)  water  depth,  and  (3)  the  depths  of  the  signal  source 
and  of  the  receiver.  In  computing  the  propagation  loss  for  the  bottom- 
reflected  path,  the  computer  program  assumes  a flat,  single-layer,  zero-loss 
(totally  reflecting)  ocean  bottom.  NISSM  II  also  provided  an  i ccurate 
estimate  of  the  grazing  angle  for  each  bottom-reflected  arrival,  as 
discussed  below. 

An  analog  reconstruction  of  the  digitally  sampled  data  (after  prepro- 
cessing filtering)  is  shown  in  the  lower  trace  of  figure  4,  where  the  four 
arrivals  resulting  from  the  detonation  of  a single  explosive  are  clearly 
separated  in  time.  The  direct,  bottom,  and  surface-reflected  arrivals  are 
isolated*  via  a UNIVAC  1108  computer  program  by  using  an  energy  threshold 
crite *ion.^ » The  program  then  computes  the  time  separation  and  pulse 
duration  of  the  various  arrivals.  Furthermore,  each  arrival  is  partitioned 
and  increased  in  duration  by  appending  zeros  to  it.  The  time  extension  of 
the  digital  data  is  required  to  satisfy  the  response  time  of  the  recursive 
digital  filters  centered  at  the  lower  frequencies.  Following  partitioning, 
the  energy  in  each  arrival  is  computed  by  squaring  and  summing  its  digital 
time  series  after  filtering  for  the  frequency  band  of  interest. 


Bandwidths  and  Center  Frequencies 

The  center  frequencies  of  the  recursive  filters  were  chosen  to  match  the 
peaks  of  the  measured  direct  path  energy  spectrum.  A measured  energy  spec- 
trum, which  was  averaged  from  20  SUS  charges  detonated  within  909  ±26  ra, 
is  shown  in  the  upper  half  of  figure  5 with  th(  center  frequencies  and  band- 
widths  of  the  recursive  bandpass  filters.  The  geometric  center  frequency, 
bandwidth,  and  upper  and  lower  effective  frequency  limits  for  each  recursive 
filter  are  listed  in  table  4.  The  digital  recursive  filters  were  five-pole 
Buttcrworth  filters. 

The  AUTOBUOYS  used  in  this  experiment  were  outfitted  with  processing 
frequency  pre-emphasis.  In  the  experiment  conducted  prior  to  MAINLOBE,  the 
AUTOBUOYS  had  been  used  for  ambient  noise  measurements  in  which  amplifier 
gain  variation  with  a frequency  range  in  the  20  to  1000  Hz  regime  is  use- 
ful. Time  did  not  permit  the  removal  of  the  frequency  pre-emphasis 

*The  upper  water  refracted  arrival  (resulting  from  a path  totally  refracted 
r^eer  the  ocean  surface)  was  also  isolated  for  identification,  but  was  not 
processed . 
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cifcuif;  prior  to  the  MAINLOBE  meatturements . The  efteot  of  the  pre-emphasis 
amplifier  has  bc'eu  removed  for  the  spectrum  shown  in  the  lower  half  of 
figure  5.  The  combined  effects  of  the  preemphasis  amplifier  and  the 
alternating  current  line  interference  filtering  resulted  in  an  effective 
broadband  of  90  to  790  lU.  Note  that  although  an  anti-aliasing  analog 
filter  with  a 750  Hit  cutoff  frequency  v as  used,  the  pre-emphasis  amplifier 
effectively  increased  the  cutoff  frequency  to  790  Hi. 


Table  4.  Filter  Parameters  Used  in  Narrowband 
Bottom  l.oss  Analysis 


(le  liter 


Frequency  (Hz) 


Bandwidth 

(oct) 


Frequency  Bounds  of 
3-dB  Passband  (Hz) 


Accuracy  of  Propagation  l.oss  and  Bottom/ 

Grazing  Angle  Determination 

As  mentioned  previously,  the  direct-path  and  bollom-ref lected  path 
propagation  losses,  together  with  the  corresponding  bottom  grazing  angle  for 
each  explosive  detonation,  were  determined  by  use  of  the  NISSM  ll  computer 
ray-tracing  program.  To  rc'duce  the  requirt'd  amount  of  computer  storage  and 
program  processing  time,  polynomial  exprt'ssions  wt're  fitted  to  the  NISSM  11 
output  tabulations.^  Polynomials  were  generated  as  a function  of 
horizontal  rang*'  between  the  signal  source  and  the  rt'ct'iver  for  (1)  the 
propagation  loss  of  each  path  and  (2)  tli*'  grazing  angle  for  tlu'  bottom 
interacted  path.  To  obtain  two  independent  determinations  of  horizontal 
range  separation,  two  additional  polynomial  functions  were  generated.  The 
first  was  a function  of  time  sepa'ation  between  the  direct  and  bottom 
interacted  arrivals  and  the  second  was  a function  of  time  separation  between 
tlie  direct  and  surfiice  reflected  arrivals. 


Thus,  for  each  detonation,  two  independent  estimates  of  range  were 
dctermliu'd . These  range  estimates  cor  i esptuided  to  tin*  t litu'  fu'parat  Ion  Inter- 
vals of  the  illrei't  path  and  surface-ref lect.eil  path  arrivals  anil  the  direct 
path  and  bottom-reflected  path  arlvals.  The  arithmetic  mean  ot  the  tV\io  range 
estimates  was  then  used  in  estimating  the  remaining  parameters  (such  as 
grazing  angle)  from  the  evaluation  of  the  corresponding  polynomial 
expression.  For  95  percent  of  the  data,  the  difference  between  tlie  mean 
horizontal  range  and  its  two  associated  values  was  less  than  ±5  percent. 

The  remaining  5 percent  of  tlie  data  had  no  more  than  110  percent  range 
error.  These  differences  result  from  differences  between  the  actual  and  the 
assumed  detonation  depth. 
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For  vungo  diftorences  of  less  than  ±5  peic«snt,  the  maximum  variation 
in  bottom  grating  angle  was  ±0.fe  deg  and  the  maximum  variation  in  propaga- 
t i on  loaa  was  tO.2  dh.  For  range  differences  of  leas  than  ±10  percent, 
tin*  maximum  variation  in  bottom  grating  angle  waa  ±1.3  deg  and  the  maximum 
variation  in  propagation  loss  was  ±0.4  dB. 

The  above  procedure  waa  modified  at  rangea  where  the  coper  water- 
refracted  arrivals  could  no  longer  be  time  i4.olaced  from  the  surface-related 
arrival.  When  this  occurred,  range  waa  computed  only  from  the  aeparation 
between  the  direct  and  bottom  pulses.  3‘his  situation  aroae  only  for  hori- 
zontal range  separations  greater  than  28  km,  corresponding  to  reflection 
grazing  angles  of  less  than  ‘).S  deg.  (lenerally,  the  upper  refracted  path 
app»>ared  only  at  horizontal  ranges  greater  than  24  km. 

3’ln'  source  It'vel  for  eai  h SU8  charge  tietonation  Wms  compnteil  by  adding 
the  direct  path  arrival  broadband  energy  to  the  value  of  energy  lost 
ali'ng  the  illi'i'ct  path  aa  prt'dlcte<l  bv  tllSSM  ll.  Tin*  predict  Ion  I rom  NISSM  1 1 In 
v\sed  in  conjunction  with  a broadband  energy  measurement;  lu'wevei’,  Ibis  does 
not  present  a prolilem  becansi*  sensitivity  studies  showed  that  there  is 
little  fiegnency  dependence  in  propagation  predictions  throughouf  the  40  to 
740  Hz  hand  over  the  range  of  0 to  3b  km  for  the  source/ reci> ivt'i  geometry 
I'mployed.  If  the  Sl'b  source  level  were  constant,  the  computed  source  level 
would  also  be  constant  with  range,  providing  that  the  propagation  prediction 
is  correct.  For  this  experiment,  the  average  stamlard  deviation  of  the 
computed  source  level  was  less  than  ±0.8  dB  for  all  of  the  data.  ll  was 
necessary  to  correct  the  NISSM  1 1 propagation  prediction  over  a brief  range 
interval  for  stations  BL-1  and  -3.  This  range  generally  occiirrt'd  around  the 
region  where  the.  direct  path  .urlvedat  the  AUTIUIOUY  at  approximately  0‘' , tl\e 
path's  tvirnlng  region.  Tt\ls  Interval  was  generally  It'ss  tlutn  2 km  amt  only 
lu'curred  at  ranges  greater  than  22  km.  The  correction  was  sutde  by  altering 
the  NISSM  ll  predlc.  Ion  iiver  the  pri'htem  range  Inti'ival  until  t lu>  I'omputed 
source  level  was  st  at  1 st  t ca  I 1 v eipial  to  mean  sourci'  levi'ls  v'omputt'd  at  tin* 
otht'i'  ranges. 


Utn'TOM  LOSS  RF.SIU.TS 


A i'omplete  si't  ot  luntom  loss  verstis  grazing  angU*  curvi's  ter  the  data 
obtained  at  all  three  statloiis  Is  contained  In  appendix  B.  Certain  aspects 
of  the  results  h.ave  been  prevhuisly  reported . For  the  purpi'se  of  cissplet 
ness  they  are  Included  In  the  ‘.ol  lowing  discussion.  To  provide  Insight 
Inti'  the  natun*  aiui  charai't  er  Ist  Ics  ef  the  bottom  loss  curves,  tin*  time 
series  structur»'  ol  the  bt>ttom  Interaction  arrival  Is  illscusseil  tlrst. 


NATLIRK  OF  THK  BOTVOM  ARRIVALS 

The  complicated  time  series  structv\re  of  u received  tixplosive  signal 
that  has  interacted  once  with  the  oceuit  bottom  can  be  simplified  through  the 
use  of  deconvolution.  The  deconvolution  process  is  used  to  determine  the 
ocean  bottom  impulse  response,  thereby,  eliminating  interference  effects  of 
the  explosive's  bubble  pulses. 
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The  method  of  frequency  domain  deconvolution  has  been  applied  to  the 
analysis  of  bottom  interacted  signals  when  the  acoustic  source  it.  ar.  explo- 
sive. 8,10,12,13  ^ simplified  mathematical  description  of  the  deconvolu- 

tion process  employed  at  NUSC  is  represented  by 


h(t)  ■■  impulse  response  of  the  ocean  bottom, 

X(f)  ■ Fourier  transform  of  input  signal  (direct  arrival), 

Y(f)  “ Fourier  transform  of  output  signal  (bottom  arrival), 

A ■ additive  white  noise  constant,  and 
G(f)  “ Gaussian  filter  function 

and  where  the  operators  F“^  { • } and  X*  denote,  respectively,  the 

inverse  Fourier  transform  and  the  complex  conjugate  of  the  Fourier  transform 
of  the  input  signal.  As  described  in  a recent  report,^  a Gaussian  filter 
function  was  used  in  the  process  to  improve  the  SNR  and  a small  fraction  of 
white  noise  was  added  to  the  direct  arrival  spectrum  (denominator)  to 
enhance  the  stability  of  the  transfer  function.®  Since  the  Fourier  trans- 
form of  the  impulse  response  represents  the  frequency  function,  the  base  ten 
logarithm  of  the  magnitude  squared  of  the  quantity  within  the  brackets  of 
equation  (2)  is,  in  actuality,  the  ocean  bottom  loss  as  a function  of  fre- 
quency. However,  to  obtain  a value  of  the  transfer  function,  this  quotient 
must  be  normalized  by  a factor  equivalent  to  the  difference  in  propagation 
losses  associated  with  the  two  signals  (i.e,  the  direct  and  bottom  arrivals). 


The  effect  of  deconvolution  can  be  seen  in  figure  6,  The  upper  trace  is 
representative  of  a typical  noise-filtered  reception.  The  shock  and  bubble 
pulses  of  the  direct  arrival  are  clearly  defined.  Because  of  bubble  pulse 
interference,  it  is  possible  to  identify  only  the  reflection  from  the  ocean 
bottom  interface  boundary.  However,  in  the  deconvolved  time  series  shown  in 
the  lower  trace,  interference  from  bubble  pulses  has  been  eliminated  by 
transforming  the  acoustic  waveform  of  the  direct  arrival  of  the  explosive 
source  into  a delta  function.  Thus,  the  direct  path  arrival  appears  as  a 
single  spike.  As  such,  additional  receptions  (representing  the  subbottom 
response  to  a delta  function  impulse  stimulus)  can  be  seen  in  this  decon- 
volved trace. 


An  approach  to  identifying  the  mechanisms  that  produce  the  subbottom 
arrivals  is  shown  in  figure  7.  Tn  this  figure,  the  time  series  of  bottom 
interacted  arrivals  are  displayed  for  eight  consecutive  values  of  the  bottom 
interface  reflection  angle.  Consecutive  receptions  are  "stacked"  along  the 
ordinate  as  a function  of  water-sediment  interface  grazing  angle.  The 
unprocessed  data  are  shown  on  the  left  and  the  deconvolved  data  are  shown  on 
the  right.  The  deconvolved  signals  were  low-pass  filtered  by  a Gaussian 
filter  with  3-dB  cutoff  frequency  at  250  Hz.  Although  the  existence  of 
subbottom  arrivals  is  evident  in  the  unprocessed  data,  their  presence  is 
distorted  by  interference  from  the  bubble  pulses.  With  the  elimination  of 
this  interference  via  deconvolution  and  the  formation  of  the  time  history 
display,  the  trend  in  the  relative  time  displacement  of  the  subbottom 
arrivals  as  a function  of  grazing  angle  becomes  more  apparent. 
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Subbottom  Reflection 

The  first  subbottom  arrival  (see  figure  7)  was  thought  to  be  a reflec- 
tion from  a subbottom  layer.  If  the  layer  thickness  is  assumed  constant  and 
independent  of  range  and  the  sound  speed  within  the  layer  is  also  assumed 
constant,  then  the  layer  thickness  and  sound  speed  of  the  layer  can  be  computed. 
Let  "At"  be  the  travel  time  difference  obtained  by  subtracting  the  travel 
time  of  the  signal  reflecting  from  the  top  of  the  layer  (the  water/sediment 
interface)  from  the  travel  time  of  the  signal  reflecting  from. the  bottom  of 
the  layer.  Let  "0g"  be  the  grazing  angle  at  the  water/sediment  inter- 
face. When  the  sediment  layer  is  thin  with  the  respect  to  the  water  layer, 
it  can  be  shown  that 


zn  1 f zn  1 Zn 

7^1  - cos  0 , 


where 


h ■ sediment  layer  thickness, 

■ sound  speed  of  water  at  the  water,/ sediment  interface, 
Cg  ■ sound  speed  of  sediment  layer,  and 
Og  ■ grazing  angle  at  water/sediment  interface. 


Letting 


y - (A  t)2, 


X ■ cos  9 , and 

g 

’ 

Equation  3 can  then  be  expressed  a:^  the  linear  equation, 

t 

y ■ mx  + b.  (/ 

Figure  8 shows  measured  values  of  (At)2  versus  cos^  0g  for  45 
points  at  site  BL-1  over  a range  of  16.9  to  47.6  deg  grazing  angle.  A least 
squares  linear  fit  was  computed  and  resulted  in  values  of  h *■  20.6  m and  Cg 
■ 1580.4  m/s,  A comparison  of  these  results  with  those  from  the  boomerang 
cores  and  archival  data  will  be  discussed  later  in  this  report. 


TR  5781 


Subbottom  Refraction  (see  figure  7) 

The  second  subbottom  arrival  behaved  in  a manner  peculiar  to  a refrac- 
ted' -ather  than  reflected,  signal.  The  refracted  signal  did  not  occur  at 
the  I'igher  grazing  angles,  but  emerged  at  26  deg  grazing  and  converged  with 
the  iirnt  subbottom  reflected  signal  as  the  grazing  angle  decreased.  This 
behavior  is  characteristic  of  an  acoustic  signal  traveling  along  a path 
through  a medium  whose  SSP  increases  with  depth. 

The  LOGO  wide  angle  seismic  measurements  made  during  the  experiment 
verified  the  existence  of  a sound-speed  gradient  in  the  sediment.  Bryan  and 
Houtz^  report  that  the  velocity  in  the  upper  250  m of  sediment  can  be 
expressed  as 


V(t)  - Vo  + Kt, 


where 


1,52  km/ s , 


K ■ 2 km/s2  , and 

where  t is  the  one  way  travel  time  of  an  acoustic  signal  propagating  down- 
ward at  normal  incidence  from  the  ocean/sediment  interface.  Equation  (5) 
can  be  reexpresscd  so  that  sound  speed  is  a function  of  depth  rather  than  of 
time: 


where  z depth, 


o 


When  the  term  — is  small,  equation  (6)  may  be  expanded  as  a 
o 

Taylor  aeries,  with  terms  of  order  2 and  higher  removed: 

V(z)  ■ Vq  + gz 

where 
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For  the  values  previously  given,  ec,uation  (7)  is  a reasonable  approximation 
(2  percent  error  in  sound  speed  to  a depth  of  300  m) . 

Travel  time  differences  between  the  refracted  and  ocean/sediment  inter- 
face reflected  signals  were  measured  for  29  points  over  a grazing  angle 
range  of  14  to  26  deg  at  station  BL-1.  These  values  were  compared  with 
predictions  obtained  using  the  Continuous  Gradient  Ray  Tracing  System  II 
(CONGRATS  II)  ray  trace  program. The  SSP  was  given  by  equation  (5), 
where  Vq  ■ 1.52  km/ s and  K ■ 2 km/s^.  Although  modeled  prediction 
agreed  reasonably  well  with  the  data,  a better  fit  was  obtained  by  account- 
ing for  the  previously  discussed  constant  velocity  layer  and  modifying  the 
parameters  Vq  and  K.  The  modified  values  of  K and  Vq  are 

K 2.16  km/s^, 

Vq  = 1.6  km/s,  and,  therefore, 

g ■ 1.35  s"^ . 

Figure  9 shows  the  measured  and  predicted  refraction/ ref  lection  travel  time 
differences  as  a function  of  reflection  grazing  angle. 


Frequency  Dependence  of  Apparent 
Acoustic  Sediment  Layers 

The  above  description  of  a relatively  simple  sediment  structure  consist- 
ing of  a 20.6  m constant-velocity  layer  overlaying  a sediment  whose  speed 
increases  linearly  with  depth  was  obtained  from  analysis  of  deconvolved 
sediment  interacted  signals  that  were  low-pass  filtered  at  250  Hz.  As  the 
frequency  passband  is  widened,  the  apparent  acoustic  sediment  structure 
becomes  more  complex.  This  observation  is  illustrated  in  figure  10. 

Filtered  deconvolved  sediment  interacted  signals  are  shown  for  five  reflec- 
tion grazing  angles.  The  signal  from  the  first  subbottom  reflector,  which 
is  clearly  seen  when  filtered  at  f^.  ■ 500  Hz,  reduces  to  background  noise 
at  a cutoff  frequency  of  200  Hz.  When  the  filtering  frequency  is  250  Hz, 
signals  from  two  subbottom  reflectors  can  be  seen.  However,  as  the  grazing 
angle  decreases,  signals  from  the  first  subbottom  reflector  become  weaker 
and  vanish  into  background  noise,  while  those  of  the  second  reflector  remain 
significant.  This  is  the  reflector  associated  with  the  20.6  m constant- 
speed  layer  previously  described. 


QUALITY  OF  DATA 

The  quality  of  experimental  measurements  can  be  degraded  by  systematic 
and  random  errors.  Systematic  error,  generally  characterized  by  biased 
measurements,  can  be  minimized  through  prudent  experiment  design  and  data 
analysis.  Random  error,  however,  results  from  the  presence  of  noise  in  the 
data  and  is  usually  unavoidable. 
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The  experimental  methodology  requiring  time-isolated  bottom  and  direct 
path  arrivals  eliminates,  or  greatly  reduces,  certain  forms  of  systematic 
and  random  errors  inherent  in  the  measurement  of  bottom  loss.  In  particu- 
lar, errors  caused  by  multipath  interference  and  improper  source  level  esti- 
mates are  eliminated. *5, 16  Errors  attributed  to  improper  estimation  of 
SUS  detonation  can  also  be  significantly  reduced. 


One  measure  of  the  magnitude  of  random  error  In  the  data  measurement 
is  the  SNR,  defined  here  as: 


SNR 


10  log 


(8) 


where 

Ig  •»  intensity  of  signal  and 
Ijq  “ intensity  of  noise. 


because  tlie  direct  and  bottom  interacted  signals  can  never  he  truly 
separated  from  noise  interference,  a modified  SNR  (SNRj,,)  is  defined  as 


SNR 

m 


10  log^Q 


+ N 


(9) 


where 

Ig+N  ■ intensity  of  signal  plus  noise 

A noise  sample  was  recorded  prior  to  the  detonation  of  each  SlIS  charge. 
From  these  samples,  Intensity  levels  were  computed  in  all  the  processing 
frequency  bands.  Thus,  for  ail  bottom  loss  measurements,  a signal-plus- 
nolse  to  noise  ratio  could  be  computed  in  all  frequency  bands  for  the  direct 
and  bottom  Interacted  arrival.  The  SNR^,  for  the  bottom  interacted  arrival 
was  always  less  than  that  for  the  direct  path  arrival.  For  the  processing 
bands  centered  at  177,  252,  325,  AA5,  600,  and  7A5  Hz,  the  bottom-interacted 
arrival  SNR„^  was,  on  the  average,  11  dB  or  greater.  The  SNR„,  was 
approximately  9 dB  in  the  broadband.  For  the  narrowband  centered  at  100  Hz, 
the  average  SNR,,j  was  less  than  A dB.  This  lower  value  results  because  (1) 
the  source  signal  spectrum  level  decreases  rapidly  for  frequencies  less  than 
the  bubble  pulse  frequency  of  151  Hz  (see  figure  5),  (2)  the  noise  spectrvim 
increases  with  decreasing  frequency,  and  (3)  there  is  a large  variance  in 
the  low  frequency  noise  resulting  from  the  short  time  duration  (approxi- 
mately A5  ms)  of  the  arrivals.  With  the  exception  of  the  100  Hz  data,  the 
relatively  good  SNRs  are  demonstrated  by  the  "closeness  of  fit"  of  the 
individual  data  points  to  the  three  points  moving  average  curve  shown  in  the 
bottom  loss  data  in  appendix  B. 


r'rtv  I,  .1. 
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GEOMETRY  DEPENDENCE  OF  BOTTOM  LOSS 

The  effect  of  source/receiver  geometry  on  bottom  loss  measurements  can  be 
observed  in  the  data  acquired  at  station  BL-1.  A receiver  (AUTOBUOY  No.  1) 
at  3638  m and  a second  receiver  (AUTOBUOY  No.  2)  at  3200  m recorded  the 
acoustic  signals  from  the  explosive  charges  detonated  at  91A  tt.  Simultaneous 
water/sediment  grazing  angle  coverage  exists  over  the  range  of  8 to  34  deg. 

In  figure  11,  bottom  loss  computed  from  the  3200  m receiver  (solid  line)  and 
the  3658  m receiver  (dashed  line)  are  compared  for  two  frequericy  bands.  In 
both  bands,  the  correlation  between  receivers  is  excellent.  In  both 
bands,  between  the  8 to  30  deg  range  the  two  curves  are  identical  to  a high 
degree  of  statistical  confidence.  This  agreement  is  significant  because  for 
each  data  point,  at  a given  grazing  angle,  th^  bottom  loss  from  the  lower 
buoy  was  acquired  using  a different  SUS  charge  from  that  used  to  acquire  the 
bottom  loss  from  the  upper  buoy.  Furthermore,  the  area  of  sediment  inter- 
action resulting  in  the  bottom  loss  data  point  at  a given  grazing  angle  is 
different  for  each  buoy.  Thus,  the  bottom/subbottom  structure  seems  to  be 
independent  of  range  for  at  least  short  distances  (4.6  km  or  less).  Also, 
the  travel  time  dlfferenoea  between  the  subbottom  and  bottom  arrivals  are 
nearly  the  same  fov  both  buoys.  An  analyses  of  ray  angles,  using  the 
model  profile  given  In  figure  9,  confirmed  that  the  bottom/subbottom 
arrival  structure  was  nearly  identical  between  the  two  buoys  for  a given 
reflection  grazing  angle.  Thus,  a geometrically  dependent  bottom  losn, 
to  the  extent  previously  reported^  (where  the  separation  between  receivers 
was  a greater  percentage  of  water  depth),  did  not  occur. 


FREQUENCY  VARIATIONS  OF  BOTTOM  LOSS 

An  example  of  typical  frequency  variations  of  bottom  loss  versus  grazing 
angle  is  shown  in  figure  12.  The  dashed  line  represents  the  3 point  moving 
average  at  a center  frequency  of  177  Hz  and  the  solid  line  represents  the  3 
point  moving  average  at  a center  frequency  of  745  Hz.  For  all  three  stations 
the  bottom  loss  generally  decreased  with  increasing  frequency  for  grazing 
angles  of  25  dog  or  greater,  and  tended  to  increase  with  increasing  frequency 
for  grazing  angles  of  25  deg  or  less. 

The  oscillations  in  bottom  loss  versus  grazing  angle  seen  in  figures  11 
and  12,  and  throughout  the  narrowband  bottom  loss  curves  in  appendix  B,  are  a 
result  of  spatially  complex  interference  patterns  created  by  the  interaction 
of  acoustic  energy  with  the  ocean  bottom.  The  broadband  processed  data 
(90-790  Hz)  tends  to  average-out  the  narrower  frequency  band  effects;  there- 
fore, the  oscillating  patterns  are  greatly  reduced. 


NEGATIVE  BOTTOM  LOSS 

The  3 point  moving  average  curve  centered  at  177  in  figure  12  shows  the 
occurrence  of  negative  bottom  loss  over  the  grazing  range  of  8 to  20  deg. 
Negative  bottom  loss  was  observed  in  about  75  percent  of  the  3 point  moving 
average  curves  computed.  As  can  be  seen  in  appendix  B,  with  the  exception  of 
data  centered  at  100  Hz,  the  level  of  negative  bottom  loss  rarely  exceeds 
-2  dB,  occurs  only  intermittently,  and  generally  is  not  observed  for  grazing 
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angles  of  20  deg  or  greater . The  observance  of  negative  bottom  loss  has  been 
reported  recently  in  other  ocean  bottom  acoustical  investigations  and  has 
been  attributed  to  the  use  ot  the  Rayleigh  reflection  model  in  the  reduction 
of  the  data.®*^  This  model  does  not  account  for  both  refracted  and  non- 
plane wave  reflected  path  signals  that  have  been  shown  to  exist  in  the 
data.®*^*^^  This  interaction  was  previously  described  and  is  illustrated 
in  figures  7 and  10.  For  the  data  reported  here,  negative  bottom  loss  is 
apparent  in  the  narrow  bands,  but  is  insignificant  in  the  broadband  where 
the  effect  has  been  averaged  out  because  of  the  frequency  diversity  of  the 
wider  band. 


GEOGRAPHICAL  VARIATIONS  OF  BOTTOM  LOSS 

A comparii'on  of  the  broadband  bottom-loss  versus  grazing  angle  data 
obtained  at  all  three  geographical  locations  is  presented  in  figure  13.  At 
station  BL- 1 , the  data  are  representative  of  the  combined  results  of  the  two 
AIJTOBUOY  receiving  systems.  The  bottom-loss  values  at  the  three  locations 
are  in  excellent  agreem(‘nt  over  the  entire  measurement  range.  Because  of 
this  close  agreement,  the  broadband  data  from  all  three  locations,  repre- 
senting 394  data  points,  were  smoothed  and  fitted.  These  data  are  shown  in 
figure  14,  where  the  solid  line  represents  a 6th  order  least  square 
polynomial  fit  to  the  394  data  points.  The  dots  represent  the  data  averaged 
in  5-deg  bins  centered  at  2.5  ♦ 5n  deg,  n “ 1,2, ...12.  The  averaging  was 
computed  in  linear  space.  Upper  and  lower  bounds  of  one  standard  deviation 
have  been  applied  to  the  averaged  data.  The  largest  value  of  the  standard 
deviation  was  l.l  dB.  Thus,  the  broadband  bottom  loss  is  characterized  by  a 
loss  of  5 dB  at  62.5  deg,  gradually  rising  to  a maximum  of  about  7 dB  at 
37.5  deg,  and  then  smoothly  decreasing  to  a minimum  of  approximately  1.3  dB 
around  7,5  deg. 

There  was  less  agreement  in  the  narrowband  bottom-loss  data  between  the 
three  locations.  This  is  attributed  to  slight  changes  in  subbottom  struc- 
ture affecting  the  complex  frequency  relationships  in  the  interaction  of 
bottom  and  svibbottom  arrivals.  The  comparison  of  bottom  loss  at  177  and  745 
Hz  for  all  three  locations  is  shown  in  figure  15.  Variations  of  up  to  10  dB 
can  be  observed. 


PRELIMINARY  MODELING  AND  INTERPRETATION 

Two  separate  methodologies  can  be  applied  to  the  problem  of  modeling  the 
bottom-loss  measurements  for  the  Hatteras  Abyssal  Plain.  Since  both  empiri- 
cal and  geophysical  models  have  advantages  and  disadvantages,  this  dual 
approach  is  necessary. 

EMPIRICAL  MODELING  VIA  POLYNOMIALS 

The  goal  of  empirical  modeling  is  to  approximate  the  data  by  a numerical 
fvinction.  This  method  can  usually  be  implemented  in  a straight  forward 
manner.  However,  it  offers  no  explanation  of  the  mechanisms  that  cause  the 
bottom  loss.  Thus,  the  model  cannot  be  extrapolated  beyond  the  range  of  the 
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data.  The  bottom  loss  data  in  this  report  have  been  fitted  to  6th  order 
polynomials.  Tabulations  and  illustrations  of  the  resultant  polynomial 
exp''essions  are  contained  in  appendix  C.  The  polynomials  allow  the  modeler 
to  compute  bottom  loss  in  decibels  as  a function  of  grazing  angle  for  each 
frequency  band  of  each  buoy  at  each  station,  for  the  combined  buoys  at 
station  BL-1,  and  for  the  combined  data  of  the  three  stations  (BL-1,  -2,  and  -3). 


INPUT  PARAMETERS  FOR  A 
GEOPHYSICAL  MODEL 

A geophysical  model  attempts  to  simulate  mechanisms  causing  bottom 
loss.  The  model  requires  a theory  adequate  to  explain  the  interaction  of 
acoustic  pressure  In  the  bottom  and  subbottom.  These  types  of  mod'^ls  require 
Input  parameters  such  as  the  sediment  thickness,  density,  attenuation,  sound 
speed,  and  sound-speed  gradient  of  each  effective  subbottom  layer. 

Some  of  the  parameters  that  could  be  used  to  model  the  bottom  loss  at 
station  BL-1  have  been  derived.  A brief  discussion  of  these  parameters  and 
variations  in  their  measurement  follows. 

The  central  region  of  the  Hattcras  Abyssal  Plain  is  characterized  by 
silts  inters trat if ied  with  clay.^®  The  total  sediment  overlaying  the 
basement  is  approximately  2.1  km,  which  is  considerably  thicker  than  the 
average  1-km  sediment  thickness  in  the  North  Atlantic. 

Earlier  in  this  report  (see  figure  9)  a partial  geophysical  description 
of  the  subbottom  was  given,  describing  it  as  consisting  of  a layer  of  con- 
stant velocity  sediment  overlaying  a layer  where  sound  speed  increased 
linearly  with  depth.  The  speed  of  sound  computed  for  the  first  layer  was 
1580  m/s.  For  the  same  area,  Horn,  et  al.,  report^®  an  average  sound 
speed  computed  from  a mean  grain  size  of  1572  m/s  (corrected  for  depth  and 
temperature)  in  the  upper  3.5  m of  sediment.  The  two  values  of  sound  speed 
averaged  from  the  boomerang  cores  to  the  north  of  station  BL-1  were  higher 
(1617  and  1609  m/s).  The  average  sediment  densities  measured  from  the  cores 
were  1.13  amd  1.30  gm/cm^.  Horn  et  al.  predicted^®  a density  of  1.77 
gm/cm^.  The  second  sediment  layer  was  characterized  by  the  linear-with- 
depth  SSP; 


V(Z)  ■ Vg  + gz 


(10) 


where 

Vg  ■ 1600  m/s, 
g “ 1.35  s" ^ , and 
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where*  z Is  the  depth  In  meters  refei  enced  to  the  boundary  with  the  first 
sediment  layer.  This  description  is  valid  to  a depth  of  approximately  145  m 
below  the  water/sedlroent  Interface,  as  this  was  the  turning  depth  of  the 
deepest  refracted  ray  measured. 

The  llnear-wlth-depth  gradient  can  be  compared  with  averaged  seismic 
measurements.  As  previously  discussed,  the  seismic  measurements  describe 
the  velocity  In  the  sediment  as  linear  with  reflection  time,  T,  referenced 
to  the  water  sediment  Interface.  The  expression  for  all  seismic  tracks  was 
given  In  equation  (5).  Following  the  approximations  discussed  In  equations 
(5)  and  (6)  a K of  2 km/s^  can  be  approximated  by  a g of  1.32s”^.  Houtz^® 
has  indicated  that  a value  of  Vq  - 1.52  km/^  may  be  unrealistically  low  for 
Abyssal  Plain, turbidities  and,  hence,  alternate  values  of  ■ 1.55  km/s  and 
k - 1.75  km/s'^  should  be  used.  This  would  lead  to  an  approximate  linear 
(with  depth)  gradient  of  1.11  s"^.  Boomerang  core  data  and  deconvolved 
records  tend  to  support  the  alternate  value  of  Vq.  Seismic  measurements  are 
made  using  methods  yielding  good  results  at  depths  deep  In  the  sediments, 
but  are  less  accurate  at  and  near  the  ocean/bottom  interface.  The  best 
resolution  available  In  the  MAINLOBE  seismic  measurements  is  approximately 
200  m.  The  paramotors  derived  from  the  deconvolved  signals  have  much  finer 
resolution  (roughly  5 m) , but  are  only  valid  to  about  145  m below  the  ocean/ 
sediment  Interface. 


SUMMARY  AND  CONCLUSIONS 

Prc  ect  MAINl.OBE  ocean-bottom-loss  experiments  conducted  by  NliSC  and 
l.DCO  In  the  Hatteras  Abyssal  Plain  of  the  Atlantic  Ocean  have  been  reported 
and  the  results  have  been  presented  in  their  entirety.  Five  significant  ob- 
servations arising  from  the  analysis  of  these  results  are  that: 

• Bottom  loss  measuieraents  were  acquired  using  a self-calibration  technique 
in  which  direct  and  bottom  arrivals  were  time  Isolated  from  all  other 
arrivals,  thereby  precluding  multipath  Interference.  As  a result,  the 
data  exhibited  excellent  agreement  over  adjacent  but  independent 
measurements . 

• Time  coincident  interaction  of  the  ocean-bottom  reflection  and  refrac- 
tions, which  are  not  ascribed  by  the  Rayleigh  model  in  the  reduction  of 
the  acoustic  data,  resulted  in  negative  values  of  bottom  loss  over  a 
limited  range  of  frequencies  and  grazing  angles. 

• A vertical  change  of  358  m In  receiver  depth  hud  little  measurable 
effect  on  the  bottom  loss  results  for  the  geometry  and  environmental 
conditions  encountered  during  these  experiments. 

• Broadband  (90-790  Hz)  bottom  loss  from  the  three  staf'  ns,  spaced 
over  an  interval  of  1 deg  of  latitude,  showed  considerable  similarity. 

A smoothed  curve  of  the  combined  stations  broadband  bottom  loss  had 

a standard  deviation  of  the  order  of  1 dB. 


• Narrowband  (1/3  octave  or  less)  bottom  loss  showed  less  similarity 
between  stations.  This  effect  is  the  result  of  small  changes  In  sedi- 
ment layer  structure  between  stations. 
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In  addition,  bottom-interacted  arrivals  at  station  BL-1  were  decon- 
volved to  obtain  the  sediment  Impulse  response,  resulting  In  the  Isolation 
and  Identification  of  reflected  and  refracted  arrivals.  A geophysical 
description  of  the  subbottom  was  custracted  from  studies  of  subbottom/ 
bottom  arrival  time  differences  as  a function  of  grazing  angle.  For 
frequencies  less  than  250  Hz,  this  description  consisted  of  a 21  m 
constant  sound  speed  layer  overlaying  a layer  where  sound  speed  Increases 
linearly  with  depth.  The  derived  sound  speed  gradient  (1.35  s"^)  of 
this  layer  was  similar  to  the  equivalent  gradient  (1.11  s“^)  obtained 
by  concurrent  Independent  seismic  studies. 

The  above  geophysical  description  can  be  used  as  an  input  to  a bottom 
loss  model.  The  output  of  such  a model  can  then  be  validated  by  comparison 
with  the  bottom  loss  results  reported  for  station  BL-1.  Because  of  the 
significant  contribution  of  the  subbottom  reflections  and  refractions,  it 
Is  recommended  that  the  bottom  loss  model  use  the  geophysical  description  of 
a layered  bottom  with  sound  speed  gradients. 
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^ AUTOBUOY  DEPLOYMENT  LOCATION 

~ SHIP  TRACK  DURING  BOTTOM  LOSS 
MEASUREMENT 

SEISMIC  PROFILING  TRACK 

O BOOMERANG  CORE  LOCATIONS 


Figure  1.  Operations  Area  for  Hatterus  Abyssal  Plain  Low 
Frequency  Bottom  Loss  Measurements 
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Figure  6.  Example  of  Signal  Enhancement  Obtained  Using  Deconvolution 
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MAINLOBE  STATION  BL-2  (28  30N  70  30W)  3353  METER  RECEIVER 


BROADBAND  (90  - 790  Hz)  3 POINT  MOVING  AVERAGE 


Figure  13.  Geographical  Variations  of  Bottom  Loss  Versus  Grazing  Angle  Data  (Broadband) 
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Figure  14.  Polynomial  Fitted  and  Averaged  Broadband  Bottom  Loss  Measured  at 
Three  Stations  in  the  Hatteras  Abyssal  Plain 
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Appendix  A 


SOUND  SPEED  VERSUS  DEPTH  PROFILES 


Results  of  the  deep  velocineter  casts  are  presented  in  figures  A-1 
through  A-5.  Separate  profiles  were  obtained  as  the  velocimeter  was  lowered 
(down  cast)  and  raised  (up  cast).  Note  that  the  velocimeter . casts  never 
reached  bottom  depths.  The  sound  speed  profile  (SSP)  was  linearly  extra- 
polated to  the  ocean  bottom  to  obtain  the  full  water  column  profile.  Table 
A-1  lists  the  relevant  parameter  for  each  profile. 


Table  A-I.  Sound  Speed  Profile  Parameters 


Figure  No. 

Station  No. 

Cast 

Actual  Depth  (m) 

A-1 

BL-1 

Down 

5600 

A- 2 

BL-1 

Up 

5600 

A- 3 

BL-2 

Down 

5600 

A-4 

BL-3 

Down 

5439 

A-5 

BL-3 

Up 

5439 

A-1 
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Appendix  B 

PLOTS  OF  BOTTOM  LOSS  VERSUS  GRAZING  ANGLE 
(ALL  STATIONS  AMD  FREQUENCIES) 


Table  B-1  provides  the  key  to  the  conditions  during  the  bottom  loss 
versus  grazing  angle  depicted  by  the  figures  in  this  appendix. 


Table  B-1.  Bottom  Loss  Versus  Grazing  Angle  Data  For  All  Stations  and  Frequencies 


Figure  Noi, 

Station  No. 

Geometry 
Source  Depth  (m) 

Receiver  Depth 
(m) 

Notea 

B-1  through  B-8 

BL-1 

91A 

3200 

- 

B-9  through  B-16 

BL-1 

91A 

3658 

- 

B-1 7 through  B-24 

BL-1 

91A 

3658 

Elxpanded  Scales 

B-25  through  B-32 

BL-1 

91A 

3658  ♦ 3200 

Combined  Receivers 

B-33  through  B-AO 

BL-2 

91A 

3353 

- 

B-Al  through  B-48 

BL-3 

91* 

3353 

- 

B-49  through  B-56 

BL-3 

91A 

3353 

Expanded  Scales 
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Appendix  C 

SIXTH-ORDER  POLYNOMIAL  REPRESENTATIONS 
OF  BOTTOM  LOSS  DATA 


The  sixth-order  polynomial  curves,  which  are  shown  in  figures  C-1 
through  C-32,  were  fitted  to  the  corresponding  three-point  moving  average 
curves  of  the  bottom-loss  data  presented  in  appendix  B,  At  station  BL-1  the 
data  from  the  combined  AUTOBUOY  receivers  have  been  used  in  forming  the 
polynomials  given  in  table  C-1.  The  polynomial  expression  for  the 
determination  of  bottom-loss  as  a function  of  grazing  angle  6 is  given  by 

Xq  + Xi0  + X20  2 + . . . +X^06,  (C-l) 

where 

BL  ■ bottom  loss  (dB), 

0 ■ bottom  grazing  angle  (deg),  and 

Xj  • constants  of  polynomial  expansion  (i  ■ 0,  1,  2,. ..6). 
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Table  C-1.  Constants  o£  Polynomial  Expression  for  the  Determination 
of  Bottom  Loss  Versus  Crazing  Angle 
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Table  C-1.  Constants  of  Polynonial  Expression  for  the  Determination 
of  Bottom  Loss  Versus  Grazing  Angle  (cont'd) 
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